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Abstract 

We consider a gas of fermions with non-zero spin at temperature T 
and chemical potential /j,. We show that if the range of the interpar- 
ticle interaction is small compared to the mean particle distance, the 
thermodynamic pressure differs to leading order from the correspond- 
ing expression for non-interacting particles by a term proportional to 
the scattering length of the interparticle interaction. This is true for 
any repulsive interaction, including hard cores. The result is uniform 
in the temperature as long as T is of the same order as the Fermi 
temperature, or smaller. 

1 Introduction and Main Results 

The physics of dilute gases at low temperature has received a lot of interest in 
the last couple of years, due to the recent experimental advances in studying 
these systems. Despite tremendous interest in the problem, rigorous results 
starting from first principles remain sparse, and often one has to rely on 
uncontrolled approximations to obtain quantitative information. This is 
true especially for dilute systems, where the interparticle interaction can 
not easily be taken into account using perturbation theory. Here, dilute 



Work partially supported by U.S. National Science Foundation grant PHY-0353181 
and by an Alfred P. Sloan Fellowship. 

© 2004 by the author. This paper may be reproduced, in its entirety, for non- 
commercial purposes. 



1 



refers to the case when the range of the interparticle interaction is small 
compared with the mean particle distance. 

The first, although admittedly not the most interesting, question to ask 
is for the ground state energy of the system. In Lieb and Yngvason 
devised a method for proving the relevant expression for dilute Bose gases. 
In this case, the energy per unit volume at density q is given by 47ra£)^, 
where a is the (s-wave) scattering length of the interaction potential, and 
a^Q <C 1, i.e., the system is dilute. (Units are chosen such that h = 1 and 
2m = 1, where m denotes the mass of the particles.) The corresponding 
expression for a two-dimensional Bose gas was later proved in '5^. 

Recently, it was possible to extend these methods and prove the corre- 
sponding result for fermions ,4.. That is, the ground state energy density of 
a dilute gas of spin q fermions is given by 



As before, a denotes the scattering length. The factor (1 — q~^) in the 
interaction energy results from the fact that only particles with different 
spin can exhibit s-wave scattering. The contribution from the interaction 
between particles of the same spin is of higher order in g. 

In this paper, we prove the analogue of at positive temperature. 
Our main result is Theorem ^ We work in the grand canonical ensemble, 
and consider the pressure of the gas at given temperature T and chemical 
potential /x. We will show that, for dilute gases, the effect of the particle 
interaction results in a contribution — 47ra^^(l — q~^) to the pressure, where 
Q is now the average density. This result holds for any temperature, as 
long as T is not much bigger than the Fermi temperature (for the non- 
interacting gas), given by Tp = (Gir"^ /q)"^^^ g^^^ (in units where = 1). 
The rational behind this formula is the following: for dilute gases, the effect 
of the interaction reduces to two-particle s-wave scattering, which can take 
place only between particles of unequal spin. This is just like in the zero- 
temperature (ground state) case. The effect of the temperature on this 
scattering process is negligible, since for T < Tp, the thermal wave length 
is of the same order (or greater) than the mean-particle distance. The aim 
of this paper is to make this intuition precise. 

We will now describe the system in detail. For simplicity, we consider 
here only the case q = 2, i.e., the spin 1/2 case. The extension to g > 2 
is straightforward. The Hilbert space under consideration is given by the 




g^^^ + Anag"^ (l - g ^) + higher order in (a^g) . (1.1) 
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fermionic Fock space for spin 1/2 particles, T = J'f{L^{Al;C^))- Here, 
Al denotes a cube of side length L. The Hamiltonian is the direct sum 
H = 0^=0 Hn, with Ho = 0, Hi = -A, and 

N 

Hn = Y,-A,+ Y1 v{xi-Xj) (1.2) 

i=l ^<i<j<N 

for N > 2. Here, A denotes the Laplacian with Dirichlet boundary con- 
ditions on Aj^. Units are chosen such that h = 1 and 2m = 1, where m 
denotes the mass of the particles. We note that both H and depend on 
L, of course, but we suppress this dependence in our notation. 

The pair potential v is assumed to be positive, radial, and of finite range 
Rq. It then has a finite and positive scattering length a. The scattering 
length may be defined as follows: if ip is the unique solution of the zero- 
energy scattering equation 

-Aip + ^vip = (1.3) 

subject to the boundary condition lim^^^^^ Lp{x) = 1, then a is given by 
o = lim|j.|_>oo — ^(x)) (see Appendix A in |H] for details). Note that 
we do not assume v to be integrable, our results also apply to the case of 
a hard core. Note also that for a pure hard-core interaction, the scattering 
length is equal to the range. 

Our main result concerns the pressure of the system at some given inverse 
temperature (3 = l/{kBT) and chemical potential fj,. It is given by 

P(/3, fi) = hm -J- In Tr^ exp {-p{H- fiN)) . (1.4) 

Here, denotes the number operator. It is well known that for systems 
with short range interactions the limit in (|1.4() exists and is independent of 
boundary conditions [HI llOj. 

We are interested in P(/3, fi) for low density, g, which is given by ^ = 
dP/dn, assuming the derivative exists. (Note that P{(3, n) is convex in /x and 
therefore the derivative exists almost everywhere. In particular, the right 
and left derivatives exist.) With low density we mean that the dimensionless 
quantity a'^g is small, i.e., that the gas is dilute. Note that there are two 
dimensionless quantities in this problem: a^g, measuring the diluteness, and 
the fugacity z = e^^. Small z corresponds to the (high temperature) limit 
of a classical gas, whereas for large z the system approaches its ground 
state. Since we are interested in the quantum aspects of the system, we will 
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consider the case when l/z is bounded. Another way of saying the same 
thing is that there are three length scales in the problem: the scattering 
length a, the mean particle distance Q~^^^, and the thermal wavelength 
We are interested in the regime where a <^ q ^ 0^^"^ 1 i-e-, where a 
is much smaller than q"^!"^ and q"^!"^ is comparable with, or much smaller 
than, 

Let Po(/3, ^) and f?o(/3, A*) be the pressure and density of a non-interacting 
gas of spin 1/2 fermions. They are given by 

^'o(/?,/i) = |(2vr)-=^ / dv In (l + zexp(-/?p2)) (1.5) 

and 

f^o(/3,/^) = ^^^ = 2(2vr)-3/ dp (1.6) 
d/i 7,83 1 + z ^ exp(/^p^) 

respectively. The factor 2 in front of the integrals takes the number of spin 
states into account. 

Our main result is the following. 

THEOREM 1. Let qq = QoiP,n) and let z = e^^" . For any a < 1/33 there 
exists a function Ca{z), uniformly hounded in z for hounded l/z, such that 

\P{(i,^l)-P^{P,^x) + 27,aQ^{P,^lf\<C^{z)aQl{aQl'y . (1.7) 

We remark that the non-uniformity of our bound for small z is not an 
artifact of our method of proof. In the limit z ^ one obtains a classical 
gas where liragQ certainly does not give the correct contribution of the inter- 
action at low density. In particular, this term depends on Planck's constant 
(which equals 1 in our units). 

Our proof below gives explicit bounds on the Ca{z) appearing in the 
statement of the theorem. Since the exponent a in the error term on the 
right side of (|1.7|) is far from optimal, however, we did not make the effort 
of writing down these bounds explicitly. Moreover, Cq(z) depends on the 
interaction potential v only through its range or, more precisely, through 
the dimensionless ratio Rq/o. Also this dependence could in principle be 
given explicitly. 

We note that go{(3, fJ-) can be replaced by the true density (of the interact- 
ing gas) in Eq. ()1.7|) . In fact, we have the following corollary of Theorem ^ 

COROLLARY 1. Let g±{P,n) = dP{(3, fj,)/dfj,^ denote the right and left 
derivatives of P{j3,^), respectively. For any a < 1/33 there exists a function 
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Caiz), uniformly bounded in z for bounded 1/z, such that 

The proof of Theorem ^ will actually not only show that the density of 
the interacting and non-interacting system are close, as claimed in (|1.8j) . but 
also the reduced one-particle density matrices (cf . Eq. (|4.61() ) . It is to expect 
that g±{P, /u) = go{f3, iJ,)—2TradQo{f3, fi)/dfi+ higher order terms in a'^go, but 
we do not have a proof of this claim. Nevertheless, we note that Theorem ^ 
implies a similar statement for other thermodynamic potentials. We state 
the following corollary for the Helmholtz free energy, but an analogous result 
holds for other thermodynamic potentials as well. 

COROLLARY 2. For g > 0, let f{P,g) = sup^[fig - -P(/?,/x)] denote 
the free energy density, and /o(/3, g) the corresponding quantity for the non- 
interacting system. For any a < 1/33 there exists a function Cq(x), uni- 
formly bounded in x for bounded 1/x, such that 

\fiP, q) - foiP, g) - 2T:ag^\ < C4Pg^/')ag^ {ag'/Y ■ (1-9) 



Theorem ^ (and Corollaries ^ and l2j can be extended in several ways, 
as will be explained now. For simplicity, we omit the proof of these asser- 
tions here. They can be proved by only small modifications of the proof of 
Theorem n given below, which is already quite lengthy itself. 

• Polarized gas: A term mS^^^ can be added to the Hamiltonian, where 
S^ot denotes the 3-component of the total spin of the particles, and m 
is a coupling parameter proportional to the magnetic field. This has 
the effect of adding a 'spin-dependent' chemical potential. Theorem ^ 
also holds in this case (with P and Pq depending on m, of course), 
if 2iTag'^ is replaced by Sirag^g^. Here, g^ and g-^ denote the density 
of spin-up and spin-down particles, respectively. They are given by 
derivatives of the pressure as ^dP/dfi ± dP/dm. 

• Higher spin: The case of higher spin can be treated in the same way. 
If q denotes the number of spin states, then the leading order contri- 
bution of the interaction energy per unit volume for a dilute gas is 
47ra(l — q~^)g'^, which reduces to (|1.7|) in the case q = 2. Also the 
polarized gas can be studied for spin higher then 1/2. 
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• Infinite range potentials: As already mentioned, the error term on the 
right side of H1.7|) depends on the interaction potential only through 
the ratio of its range to its scattering length, Rq/u. By cutting off 
the potential v in an appropriate (^j-dependent) way, it is therefore 
possible to extend Theorem ^ to infinite range potentials (with finite 
scattering length), with possibly a worse error term than the one given 
in (|1.7)) . (Compare with the corresponding discussion for the Bose gas 
in the appendices of ^ and |S|.) 

• The two-dimensional gas: A corresponding result can also be derived 
for a Fermi gas in two dimensions. The leading contribution of the in- 
teraction energy for a dilute gas is then 27rf3^/| In a'^g\ per unit volume. 
This was shown in for the ground state, i.e., at zero temperature, 
and the methods developed in this paper can be used to obtain this 
result also at positive temperature. 

Before giving the full proof of Theorem ^ we start with a short outline 
to guide the reader. In the next Section |2 we state some preliminaries 
that will be useful for our proofs. We introduce the pressure functional 
which defines a variational principle for the pressure. We also state some 
useful properties of the non-interacting system. Sections El and contain 
the proof of Theorem ^ This proof is divided into two parts, the lower and 
upper bounds to the pressure. Finally, in Section [5J we give the proof of 
Corollaries ^ and El 

For the lower bound to the pressure, given in Section |31 it is necessary 
to construct an appropriate trial density matrix for the pressure functional. 
As in the zero temperature calculation in j^, we find it necessary to choose 
a trial density matrix that confines particles into small boxes in order to 
control the average particle number in each box. The construction in each 
small box is done in Subsect. 13.21 We then proceed with the calculation of 
the variational pressure in Subsect. 13.31 We use similar methods as in |4j to 
estimate the energy of the state. In addition, it is necessary to estimate its 
entropy, which is done with the aid of Lemma El The price one has to pay 
for using the box method are finite size corrections, which are estimated in 
Subsect. 13.41 The final result is then stated in Subsect. 13.51 

The upper bound to the pressure, given in Section has two main 
ingredients. First, an operator inequality proved in j3i Lemma 4], which 
allows for the replacement of the interaction potential f by a "soft" potential 
U, at the expense of the high-momentum part of the kinetic energy. See 
Subsect. 14.11 Secondly, an improved version of subadditivity of entropy in 



6 



Subsect. 14.21 which was recently derived in As shown in Subsect. 14.31 
this estimate allows to prove that the reduced one-particle density matrix 
of the spin-up particles, for fixed positions of the spin-down particles, is 
close to the Fermi-Dirac distribution for non-interacting particles. This 
property of the reduced one-particle density matrix is then used to show, in 
Subsect. 14.41 that first order perturbation theory with the soft potential U 
gives the correct answer for dilute gases. 



2 Preliminaries 



Since the Hamiltonian H does not depend on the spin variables and, in 
particular, commutes with the operators counting the number of spin-up 
and spin-down particles, the problem can be reformulated in terms of two 
species of spinless fermions. More precisely, JT = J^p(TCi) T'-piJii) = 
J^^ ® J^i, where Hi = L'^{Al;C) denotes the one-particle space for spinless 
particles. We label particle coordinates in the first factor by x^ and in the 
second by x^. The Hamiltonian in this representation can then be written 
as H = ^N,M=o with 

N M N M 

i=l k=l i=l k=l 

+ Yl Y (2-1) 

l<j<j<Af l<k<l<M 

The first two terms are simply the kinetic energies of the spin-up and spin- 
down particles, and the interaction potential is divided into three parts, 
corresponding to interaction between particles of the same and of differ- 
ent spin, respectively. In a sector of fixed particle numbers N and M, 
we denote the particle coordinates collectively by = (x|,...,x]y) and 
= (xj, . . . , x^^). 



2.1 The Pressure Functional 

The pressure (|1.4|) can be computed via a variational principle. For T a 
density matrix, i.e., a positive trace class operator on with Trj^-F = 1, we 
define the pressure functional [T] by 



S[T] 



(2.2) 



7 



where S[T] = — Trjr(rinr) denotes the (von Neumann) entropy. (The ex- 
pression ()2.2() is well defined if the eigenfunctions of T are in the quadratic 
form domain of H and N. Otherwise, we can take it to be +00.) Let 
P^{P, /x) denote the maximum of 'P^[r] over all density matrices. The max- 
imum is uniquely attained by the grand-canonical Gibbs density matrix, 
given by exp{—(3{H — fiN)) /Tryr exp{— (3{H — fiN)). Hence the pressure 

P^{(3,ij) is given by P^{/3,fi) = L-^^'MnTr^ exp (^-/3(/7 - ^iV)) , and 

P(/3,/x) =limL^ooi^^(/?,/x). 

At zero temperature, i.e., when /? = 00, this variational principle reduces 
to the usual variational principle for the ground state energy. Note, however, 
that at positive temperature the functional ()2.2() is not linear in the density 
matrix. 

2.2 The Ideal Fermi Gas 

For later use, we also define the pressure functional for the noninteracting 
gas, 7^(f [r]. It is defined in the same way as 'P^[r] above, with H replaced 
by the non-interacting Hamiltonian H^^^ = 0^jv/=o -^.WM' where H^^^j^ is 
given as in 1)2.1(1 but with u = 0. We denote the (finite volume) pressure for 
the non-interacting gas by Pj" (/3, /i). It is given by 

^o'lA/") = |;^T>Hi ln(l + zexp(/3A)), (2.3) 

which reduces to (|1.5|) in the thermodynamic limit. We note that this ex- 
pression can also be obtained from a variational principle for the reduced 
one-particle density matrix (see, e.g., 12 ). Namely, 



-iL3Po^(/3,Ai) =inf 



Tr„,(-A-/x)7--^5[7] 



(2.4) 



where the infimum is over all positive trace class operators 7 on Tii with 
< 7 < 1, and S[^] is given by 

%]=Tr>^, [-7ln7-(l-7)ln(l-7)] • (2.5) 

The reason for the factor 1/2 in front of Pq are the 2 different spin states, 
which we have not accounted for in the functional. The infimum in (|2.4|1 is 
uniquely attained at 70 = (1 + z~^e~^^)~^ . 



8 



3 Lower Bound to the Pressure 



We start the proof of Theorem ^ by deriving a lower bound to the pres- 
sure. Since -P(/3, ^) is determined by maximizing the pressure functional, a 
lower bound can be derived using an appropriate trial density matrix in the 
pressure functional H2.2|) . 

3.1 The Box Method 

It will be convenient to divide space into small boxes of side length I and 
confine the particles to these boxes. By choosing I appropriately, we can 
then control the average particle number in every box. Moreover, if we keep 
these boxes separated by a distance i?o, there is no interaction between 
particles in different boxes. 

More precisely, pick an integer / and divide the interval [0, L] into / 
intervals of equal length. We choose / such that £ = L/I — Rq >0. From 
the variational principle defined by (|2.2|) we can infer that 

L3p^(/3,/i)>/W(/3,/i), (3.1) 

where the factor is the number of boxes. Dividing by Lp' and letting 
L — > oo and / ^ oo in such a way that L/I converges to some number 
greater than i?o, we see that 

for any £ > 0. 

3.2 Construction of the Trial Density Matrix 

We now construct a trial density matrix for 'P^[r]. For fixed /3 and /u, let 
^'o = ^'o(/?5/i)) and let K > Let Q be the projector onto the subspace 
of Hi = L2(A^;C) where -A < KqI'^, i.e., Q = e{Kgl^^ + A). Here, 9 
denotes the Heaviside step function, given by 

On Tq = T{QTi.i) T{QTi.i), let Hq denote the second quantization of 
—AQ, and let Tq be the corresponding grand canonical Gibbs density ma- 
trix, defined as Tq = exp{—P{HQ — fj.N))/TTj^g ex.p{— P{Hq — nN)) . We use 
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the same symbol for the density matrix on T ^ being Tg on the subspace Tq 
and on the orthogonal complement. Let qq = £~'^TrjrNTQ be the average 
density of Tq. By explicit computation, 

We can decompose Tq as 



'^Q = '^^a\'4'a){lpa\, (3.5) 
a 

with Aq, > 0, Y^a^ct = 1, and {tpa} an orthonormal set in J^q. Moreover, 
we can always choose the tpa to be products of Slater determinants of Na 
j-particles and Ma j-particles, respectively, for some Na,Ma € N. Given 
such a V'a) we define 

with F^'^^ given as follows. Pick some s > 2Rq and let g : M'^ R be 
function with < (7 < 1, having the property that g{x) = for |x| < s and 
g{x) = 1 for \x\ > 2s. We may also assume that \Vg\ < const. for some 
constant independent of s. Moreover, for some > R > Rq, let / : i— > M 
be given by f{x) = ip{x)/{l — a/R) for < R and 1 otherwise. Here, if 
denotes the solution to the zero-energy scattering equation (|1.3() . Note that 
/ is a continuous function, since (p{x) = 1 — a/\x\ for |x| > Rq. We define 



N M 

^^'^(xT,xi)= n 9ixi-x]) n 5(4--i)nn/(^' 

l<i<j<N l<k<l<M i=lk=l 



X, - X^) 



(3.7) 

As a trial density matrix for 'P^[r] we choose 

r = ^A„|./.„)(</.„|, (3.8) 

a 

with (f)a defined by (|3.6|) and (|3.7j) . Note that the are not orthogonal, 
but they are normalized, and hence Tij^T = 1. 

3.3 Calculation of the Variational Pressure 

We now derive a lower bound on the variational pressure [F] , with F given 
in (|3.8)) . We start with the expectation value of the energy. 
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which we have to bound from above. This expression can be estimated using 
the same methods as in |^ Sect. IV]. More precisely, the calculation in |^ 
shows the following. 

Lemma 1. For N,M > 0, let Di and D2 denote Slater determinants of 
N and M orthonormal eigenf unctions of the Dirichlet Laplacian on a cube 
of side length I, respectively. Let k denote the maximal kinetic energy of 
these N + M functions. Let il}{X\X^) = Di{X^)D2{X^)F^^'^'\X\X^), 
with F^'^^ given in \3. Then, for some constant c > 0, 

M||W < ,10,0,) + s™^ (1 + cN-^r^ + eM-/3 

NM s^/^a^/^ 
+ca—^E{R, s,N + M, k, i) + c{N + M)^/^ 



£4 ' 

(3.10) 



where E is the function 



E{R, s, n,k,e) = ^ + s^k + ^+ n^l^{s/lf . (3.11) 

This lemma was proved in [4^ for the special case when Di and D2 are 
Slater determinants of the lowest N and M eigenfunctions of the Dirichlet 
Laplacian, respectively (compare with |^ Eq. (46)]). The proof, however, 
goes through without change in the general case, the only difference be- 
ing the maximal value of the kinetic energy, k, which enters the bound 
through the estimate in Lemma 2]. More precisely, the value r?!"^ in 
[U Lemma 2] has to be replaced by k in the general case considered here. 

Note that the error term in Lemma ^ is not uniform in the particle 
number. For this reason, it is necessary to confine the particles into small 
boxes, as done here. 

In the case of interest in ()3.9() . the number of particles is at most N+M < 
2Tr^,Q < ce^K^/'^QofoT some constant c > 0, and the kinetic energy of each 
factor in the Slater determinants is at most k < Kg^^^. Hence Lemma ^ 
implies the upper bound 

P 



+ca'-^E{R,s,K'/'goi',KeT,i) 
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for some constant c > 0. We now insert this estimate into p.9|) . We have 
^A«iV„M, = Tt^N'^N^Tq = ill^QQ^ , (3.13) 

a 

where denotes the number operator on J^-^. Moreover, using convexity 
of X I— > x^^'^, it follows from Jensen's inequality that 

J2^aN'J'M^<{le^eQf\ (3.14) 

a 

and likewise with Na and Ma interchanged. Also J2a '^ai^a + Ma) = 
Tt^NTq = QQ^. Thus we obtain the upper bound 

+E{R, s, K^/^go£\ Kgl^^i)]) (3.15) 



for some constant c > 0. 

The next step is to calculate the average particle number. By construc- 
tion of r, 

Tr^iVr = J]A„(0«|iV|0„) = ^A«(Va|/^|Va) = Tr^^rQ. (3.16) 

a a 

It remains to derive a lower bound on the entropy of F. For this purpose, 
we need the following Lemma. 

Lemma 2. Let T be a density matrix on some Hilbert space, with eigenvalues 
Xa > 0. For {Pa} (not necessarily orthogonal) one- dimensional projections, 
let r = XaPa- Then 

S[f]>S[T]-ln\\j:^Pa\\ . (3.17) 
Proof. Using twice concavity of the logarithm, 

5[f]-5[r] = -^A,TrP„ln(A-if) 

a 

> -^A,lnTrP„A-ir 

a 

> -lnTr(s„P„f) >-ln|E„P„|| . (3.18) 
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Let X = maxa ||F^"'*^"^a||"2_ Then, 



(3.19) 



Note that F^"'*'^« depends on a only through the particle numbers and 
Ma- Denoting the sum over a sector of fixed and Mq, by and using 
the fact that the ipa are orthonormal, we see that 



This hold in every sector, and hence 

<X- 

a 

Lemma 121 thus implies that 

5[r]>5[rQ]-lnx, 



(3.20) 



(3.21) 



(3.22) 



and it remains to derive an upper bound x- Equivalently, we need a lower 
bound on the norm ||F^"'*^°''0q,||. This can be obtained as follows. 

Lemma 3. Under the same assumptions as in Lemma^ 







niin{N,M} 


> 




+ 



l-c{N + Mf/^{s/if (3.23) 



for some constant c > 0. Here, = max{t, 0} denotes the positive part. 

Proof. We write F^'^{X\X^) = Gn{X^)Gm{X^)H{X\X^), where Gtv, 
Gm and H denote the three different factors in ()3.7I) . From ^ Lemmas 1 
and 3] we can infer that 



dXUX^ Di{X^fD2{X^fGN{X'^fGM{X^?H{X\X^f 
> j dXUx^ Di{X^fD2{X^fGN{X'^fH{X\X^f 
l-cM^I'^\\A'-^^f{s/lf 



dX^ Dl{X^ fGN{X^ f{dQiKx^) \l - cM^/^\\K-^^f{s/(f 



(3.24) 
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Here, A^t denotes the M x M matrix 

("^^0 nm 



N 



(3.25) 



with ifn denoting the M eigenfunctions of the Dirichlet Laplacian that con- 
stitute the Slater determinant D2, and || • || stands for the operator norm. 
Note that because of the factor Gm the integrand in 1)3.24(1 is only non-zero 



if |x- — Xj \ > s for all i ^ j- In this case, Lemma 2 in jl] implies that 



II- A 



XT I 



< c 



(3.26) 



(Again, as already mentioned after Lemma^ the factor ri^^^/l"^ in the state- 
ment of 01 Lemma 2] has to be replaced by k in the general case considered 
here.) In particular, since < A < I, this estimate implies that 



||A 



-1 1 

XT I 



< 



and that 



detAj((^T > 



+ s^k 



c I + s^k 



M 



(3.27) 



(3.28) 



By inserting these two bounds into ((3.24|) and using again Lemma 3] to 
get rid of the Gn in the integrand this implies (|3.23() in the case M < N. 
The case N > M follows in the same way, interchanging the estimates for 
the X'^ and X-^-particles. ■ 

To apply this lemma, we use again the fact that in the case of interest 
the number of particles is at most N -\- M < 2Ttt-i^Q < c£^K^^'^qo for 

2 /3 

some constant c > 0, and the kinetic energy k is bounded by Kqq . Hence 
Lemma 01 implies that 



1 

- > 

X 



1 



aR^ 



+ s^Kq, 



2/3 



1 - (s^Qo 



5/3' 



(3.29) 



for some constant c > 0. In combination, (|3.15|) . H3.16|l and (|3.22|) imply the 
lower bound 



1 



V'[T] > V',[TQ]-27rag'Q-—lnx 



-caQQ 



I R-S.^/S -1 3« -1/2 



(3.30) 
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with X bounded by H3.29() . 



3.4 Approximating Traces by Integrals 

The pressure of Tq is easy to compute: 

V'^[Tq\ = ^Trn, In (l + Qexp ( - /?(-A - ^))) . (3.31) 

We have to compare this quantity with the true pressure of the non-inter- 
acting gas in the thermodynamic hmit, 



Po{P,^,) = ^{2^^y'^ [ dpln(l + exp(-/3(p2_^))) . (3.32) 



To this end, we note the following: 

Lemma 4. Let f : M_|_ be a monotone decreasing function, and let A 

be the Dirichlet Laplacian on a cube of side length I. Then 



(27r)-3 / dpfip^) > r^Tvn, /(-A) > (27r)-3 / dpfip^ 
Jr3 J^3 



1- 



(3.33) 



Proof. Note that the spectrum of —A is given by [(7r/^)N]'^. Considering 
the trace as a lower Riemann sum to the integral, we immediately obtain 
the first inequality. To obtain the second, we consider the trace as an upper 
Riemann sum to the integral over the region where Pi > Tr/£, 1 < i < 3, 
with Pi denoting the components of p. Hence 

3 

{2^)H-^TthJ{-^)> I dpf{p')-Y,8[ dpfip'). (3.34) 

Jr^ Jo<pi<TT/e 

Since / is monotone decreasing, we can estimate 

/ dpfip') [ dpfip') = f dp^f{p') , (3.35) 

J0<Pi<7r/e 4 £ 7]g2 8 £ J ^3 \p\ 

where the integral in the second term is over the plane M? . Combining ()3.34|) 
and (|3.35|1 we arrive at the second inequality in (|3.33j) . ■ 
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From this lemma, we immediately see that qq < qq. Moreover, 



T'oITq] - Po{P,fM) > -|(2^)-' / ,,3 dp In (1 + zexp(-/3p2)) 



(3.36) 



In the integrand in the first integral, we estimate ln(l + x) < x as well as 
\p\ < (2/?)-i/2exp(i/5p2) and obtain 

|(2vr)-=' ^^^dpln(l + zexp(-/3p2)) < _!__£_ exp (-^PKe'^' 



The second term on the right side of (|3.36j) is bounded from below by 
3 1 r . I . ,^ , n.. 3 1 



~~r^lyn / dp— r In (1 + zexp{-p'^)) = -—^—r-g{i 
A-K^ [3^1 \p\ ^ ^\ ^ jj i-jT^ p^i^y 

and therefore 



(3.37) 

r 

(3.38) 



7r2/55/2 



(3.39) 



By combining (|3.39|) with (|3.3fl|) and the estimate qq < go, we thus obtain, 
for some constant c > 0, 

P'{P,I^) > Po{f3,fi)-27ragl--^lnx 



-caQQ 



3/3 



1/2 



9{z) + A= exp (^-^PKqI^^ 



V2 



(3.40) 



^2^5/2 

with X bounded by 1)3. 29(1 . 
3.5 Final Result 

We are still free to choose £, R, s and K. We choose 

R = a{a^eo)-'^'' , s = aia^eor''/'' , £ = eo'^'i^^ 60)-''^'' (3.41) 
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and K = (a'^Qo) ^^^"^ for some e > 0. With this choice the first term in 
square brackets in (|3.4U|) is bounded by 

(3.42) 

for some c > and a? Qq small. Moreover, 

In X < const. -^473 {^qI'^) ^/^^"^^^^ . (3.43) 

2 /3 

Note that /3^q is a monotone increasing function of z and, in particular, 
l/(/9^Q ) is bounded for bounded 1/z. 

The first term in the last line of (|3.4()|) is 

^ ( \ 2 9{z) , 1/3x1/27 

-^^g{z)=aQ^ ^^^^l^^^ {aQ^ ) . (3.44) 

Now fiq^ ~ ln(z) for large z, and also g(^z) ~ ln(z) for large z. Hence the 
fraction in 1)3. 44() is uniformly bounded in z for bounded Xj z. The remaining 
term in 1)3. 4U() is 



5/3 1 r / , o., , 



2/3n5/2 



exp (^-i/^^o^Vi^T^Mnz) . (3.45) 



For large the term in square brackets is strictly less than one, again 
uniformly in z. Hence we see that the expression ()3.45() is exponentially 

small for small c? and, in particular, bounded by const, {p^ for 
any exponent p, with a constant that depends on p (and e), of course, but 
not on z (for bounded l/z). 

To summarize, we have thus shown that with the choice of parameters 
as above, (|3.4fl|) gives, for any e > 0, 

P\(i,ix) > PoiP,fi)-27ragl (l + C,(z)(a£>J/=^) , (3.46) 

with some constant C£{z) that is uniformly bounded in z for bounded 1/z. 
To complete the estimate, we have to insert this bound into (|3.2j) . Thus we 
still have to estimate 

^P,i(3,,) = ?^ael[ael'Y''\Pi'r''"m,M . (3.47) 
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Now Poil,f3n) ~ (lnz)5/2 ^j. ^^^^.g^ therefore {Pgl'^^^/'^Poil, Pfi) is 

uniformly bounded for bounded 1/z. Altogether, this implies that, for some 
constant Ce{z) uniformly bounded for bounded 

P{(3,^Ji) > Po(/3,/i) - 2^aQl (l + Ce{z){aQYY^''') ■ (3-48) 
This finishes the proof of the lower bound. 

4 Upper Bound to the Pressure 

In maximizing the pressure functional we can restrict ourselves to density 

T 

matrices that do not mix particle numbers. More precisely, if Qjy denotes the 
projection onto the sector of particles in then 'P^[r] < 'P^[r], with 

r = X^Af M ^jv^i/^^TV^M- This follows from the fact that the entropy 
is non-decreasing under this transformation |121 2.1,11.4]. Hence we can 
assume that F = ^^ym '^n,m-, with Tjsj^m (not normalized) density matrices 
for N |-particles and M |-particles, respectively. For simplicity, we may also 
assume that F is symmetric with respect to exchange of | and [. This is 
certainly no restriction in the case considered here, and leads to simpler 
expressions by shortening some of the formulas. 

4.1 Lower Bound to the Hamiltonian 

We start with a lemma which is essentially a generalization of a result by 
Dyson to bound the hard potential v from below by a soft potential 
f7, at the expense of some kinetic energy. In the following, / denotes 
the Fourier transform of a function /. We use the convention f{p) = 
(27r)-3/2/(ix/(x)e-^P^. 

Lemma 5. Let x be a radial function, < x ^ 1? with h = 1 — x ^ 
L^{R^) n L~(M3). For R > Rq, let 

fnix) = sup \h{x -y) - h{x)\ , (4.1) 

\y\<R 

and 

WRix) = ^fnix) dyfniy). (4.2) 

Let U be a positive, radial function, supported in the annulus Rq ^ \x\ ^ R, 
with Jj[j3 dx U{x) = 4tt. If yi, . . . , um denotes a set of M points in R^, with 
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\yk — yi\ > 2i? for all k ^ I, then, for any e > 0, 

M M 

-Vx{p?y + i ^ 7;(x - yfc) > ^ ((1 - £)a;7(x - yfc) - ^wr{x - yu)) 

k=l k=l 

(4.3) 

in the sense of quadratic forms. Here, x{p) stands for the multiplication 
operator in momentum space. This operator inequality holds for all functions 
in if-'^(M^) and therefore, in particular, for functions supported in the cube 
Al. 

The proof of this lemma can be found in [Ij Lemma 4 and Cor. 1]. Note 
that, by construction, either wr G L^(M^) n L°°(M^), or wji = oo identically. 
Lemma El implies the following lower bound to the Hamiltonians (|2.1|) : 

N 
i=l 

M 



+ E [-^i (i - ^(^i)') + (4)] • (4.4) 



k=l 

Here, Wy is the potential 

a 

e 



Wy{x)= [{l-e)aU{x-yj)--WR{x-yj)) , (4.5) 

{k:yk(^YR} 



where Y = (yi, . . . , yM) is any set of coordinates in and Yr C 1" is the 
subset of y^-'s whose distance to the nearest neighbor in Y is at least 2R. 
We neglect the interaction with yj's that are not in the set Yr, which can 
only lower the energy. Note that also the interaction terms among particles 
of equal spin are dropped for a lower bound. 

The expectation value of the potentials Wj^t and Wxi can be written 
in the following convenient way. For T = ®j^m^n,m a density matrix 
on (g) Ti, and X'^ = (xj, . . . ,xjyj) some fixed coordinates of the spin- 
down particles, let the operator be defined through the integral kernel 
V-^^ {X\Y'^) = Tj^^m{X\X^,Y\X^). Since Lat^m is a trace class operator, 
this expression is well defined for almost every X^ by an eigenfunction expan- 
sion of rAT^A/- The same is true for n{X^), given by n{X^) = Trj^i [®7v '^'n^\- 
Note that n{X^) is the probability density for having exactly M spin-down 
particles at the positions XK In case n {X^) > 0, let 



rf =n(xi)-i0rf . (4.6) 

N 
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This defines ri*"^ only if n{X^) is non-zero; only in this case it will be used 

below, however. Note that r^"^ is a density matrix on J^^, which can be 
interpreted as the density matrix of the t-particles for a fixed configuration 



of the |-particles. If denotes the reduced one-particle density matrix of 



r^''''^ , then the expectation value of the potential Wxi in the state T can be 
written as 



(4.7) 



I 

M- 



where we introduced the short hand notation ^ dX^ = / ' ' ' dx 
Note that f dX^ n{X^) = 1, and f dX^ n{X^) Vf^ = ee Tr^^ T. 

Under the assumption that the density matrix is symmetric with respect 
to exchange of | and |-particles, the lower bound to the Hamiltonian in 
(|4.4() can thus be written as follows: 

T^rHT>24dX^n{X^)TvnA[-V{l-x{p?W + Wxi)lf] • (4.8) 



Note that the kinetic energy term does not depend on X^ and, therefore, 
the integration affects only . Note also that ^ dX^ n{X^) '^^^ = 7^, the 
reduced one-particle density matrix for the |-particles. 

For reasons which will become clear later, we find it convenient not to 
use up all the kinetic energy in the bound (|4.8j) . however. More precisely, 
we pick some < 5 < 1 and < x < 1 and write —A as 



with 



-5 A - (1 - ,5) (1 - x)Vx(p)^V + , 



.(l_5)V(l-(l-x)x(p)')V. 



(4.9) 



(4.10) 



Applying the above estimate only to the second term on the right side of 
(|4.9|) and using positivity of the interaction potential v, we obtain that 



Tr^OT > 2Tr'H, [{-SA + h^)j^] 

+2(1 - 6){1 - >c)^dX^ n{X^)Trn, 



(4.11) 



Eq. ()4.11|1 is the final result of this subsection. To estimate this expression, 
we will show that, for any fixed X^, the one-particle density matrix 
is close to the corresponding expression for non-interacting particles (which 
does not depend on X^, of course). We do this in Subsection 14.31 below. 
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It remains to choose x- Let 77 : ^ M_(. be a smooth radial function 
with r/(x) = for |x| < 1, ri{x) = 1 for > 2, and < r}{x) < 1 in-between. 
For some s > we choose 

X{p) = r]{sp) ■ (4.12) 

The potential Wy then depends on e, a, R and s. Note that with this choice 
of X the corresponding /i = 1 — x is a smooth function of rapid decay and 
hence, by simple scaling, the corresponding potential wji defined in 1)4. 2|) 
satisfies, for R < const, s, 

-ft It 

h^R 00 < const. and h < const. (4.13) 

for some constants depending only on rj, which is fixed once and for all. 
Moreover, if ly^ — y^l > 2R for all A; 7^ /, then 

M ^ 

^ wr{x - Vk) < const. (4.14) 
fc=i 

independently of x and M. 

We are also still free to choose the potential U in Lemma El We choose 

^(^) = / 3 - i?3)-i for i?o < \x\ < R, 
\ otherwise . 

We then have the estimate 

||t^y||oo <max|-^^3^^, const.-^l , (4.16) 



independently of Y. 



4.2 Improved Subadditivity of Entropy 

For r a density matrix on J^-^- J^i, let Fj = TrjFj^ F and F| = Trjr,^ F be the 
density matrices of the subsystems of j and J,-particles, respectively. It is 
well known that the entropy 5[F] is subadditive (see, e.g., ^1 Ineq. (2.2,13)], 
i.e., 

5[F]<5[Ft]+5[Fj, (4.17) 

where the entropies on the right side are defined by taking the trace only over 
J^i and respectively. Moreover, if 7| denotes the reduced one-particle 
density matrix of Fj , then 

5[F|]<5[7i], (4.18) 
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where S[y] = Tr^^j (— 7ln7 — (1 — 7) ln(l — 7)) is given as in (|2.5|) [T^ 
Ineq. (2.5,18.5)]. Note that both KTT^ and (IHHj) are equahties if T is 
the grand-canonical Gibbs density matrix of a non-interacting system. 

We are going to need the following refinement of subadditivity of entropy. 
Its proof is given in Cor. 4]. 

Lemma 6. Let T = '^n,m be a density matrix on T-^® T^. With the 

notation introduced in Subsect. \4 1[ 

s[r] < s[r^] + ^ dx^ nix^) s[rf] . (4.i9) 

Note that f dX^ n(X^) = 1, andfdX^ n{X'')Tf^ = T^. Hence, by con- 
cavity of >S'[r|], inequality (|4.19|) is stronger than the usual subadditivity of 
entropy in 1)4. 17(1 . The last term on the right side of ()4.19|) is the average 
entropy of the j-particles for fixed |-particles, whereas S\ri^] is the entropy 
of the state of j-particles averaged over all configurations of the |-particles. 

4.3 A Priori Bounds on the One-Particle Density Matrix 

In this subsection we will show that the one-particle density matrices for 
fixed |-particles, 7^''^, are close to the corresponding expression for non- 
interacting particles, provided the state T that defines them has a variational 
pressure P^[r] close to the true pressure P^{(3,^) of the system, i.e., it is 
an approximate maximizer of the pressure functional 1)2.2(1 in a sense to be 
made precise below. We call such a bound an a priori bound. 

This subsection is split into four parts. In Part 1, we will derive a 
bound of the sort needed. This bound will not be uniform in the fugacity 
z, however, and will be useless for very large z (when the system is close 
to its ground state). For large z, however, we will use a different method 
to obtain a similar bound, by comparing the state with the ground state 
of the non-interacting system. We do this in Part 3. Before that, we use 
the method of Part 1 to compare the non-interacting system with differ- 
ent boundary conditions (Part 2). Some calculations are easier to do with 
periodic boundary conditions than with Dirichlet, hence the usefulness of 
this estimate. Finally, we give a summary of the result of this subsection in 
Part 4. 
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4.3.1 General a priori bound 

Using the fact that f > 0, we can mfer from Lemma IHl and (|4.18jl that 

- L^V'^IT] > Ti-n, [(-A - ^) 7i] - ^Shi] 

+^dX^n{X^) {rvn, [(-A-/.)7r] --^5[7f ]) (4.20) 

The term in the first hne on the right side of (|4.2n|) is bounded from below 
by -\L^P^{[j,ii) because of JTH). The same is true for the term in the 
last line, but we will need a refinement of this inequality, given in Lemma |7| 
below. 

For r an approximate maximizer of P^[r], we have an upper bound 
on the left side of H4.2U() . derived in the previous section, and therefore this 
yields an upper bound on the last expression on the right side of (|4.20|) . This 
bound can be used to get information on the one-particle density matrices 

. We need the following lemma. 

Lemma 7. Let h he a self- adjoint operator on a Hilbert space 7i, such that 
e~^ is trace class. For 7 a fermionic one-particle density matrix (i.e., a 
trace class operator on TC with < 7 < 1), define the functional 

£h[l]=TThj-S[^]. (4.21) 

Let 7/i = (1 + e^)^^ be its minimizer, and f{h) = £h[lh\ = — Tr ln(l + e~^). 
Then, for any 7, 

^^hb] >/W + 2Tr(7-70^ (4.22) 



and also 



hlJi-JK J 2|T¥(7-70|+TV7;, ^ ' 



Note that this lemma implies, in particular, that 7 —> 7/1 in trace class 
norm, if Shi'j] fih). Eq. (|4.22j) implies convergence in Hilbert-Schmidt 
norm, but because of 1)4.23(1 also the traces converge, and therefore the 
convergence is in trace class norm (see |13| lllj or Ineq. (|4.35j) below). 

Proof. We write 

£hb]-f{h)=TTg{j,jh), (4.24) 

with 

9(1, Ih) = 7ln7 - 7ln7fe + (1 - 7) ln(l - 7) - (1 - 7) ln(l - 7,^) . (4.25) 
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The function g has the integral representation (for < x,y < 1) 

g{x,y) = 



az [x — zj I — h 



z 1 — z 



(4.26) 



Using that 1/z + 1/(1 — z) > 4 in the integrand, we obtain the lower 
bound g{x,y) > 2{x — y)'^. Hence, by Klein's inequality ^1 Ineq. (2.1,7.5)], 
Trg'(7,7/j) > 2Tr(7 — 7/1)^, and (|4.22j) follows. Moreover, estimating 1/z + 
1/(1 — z) > 1/z > 1/ max{x, y} > l/{\x — y\+y) in the integrand in 1)4. 26() . 
we obtain 



gix,y) > 



1 {x-yy 



2 sup [b{l-b)\x-y\ -b^y] . 



(4.27) 



2\x-y\+y o<fe<i 

Hence, again by Klein's inequality, 

Tr5(7,7/.) > 26(1 - 6)|TV(7 - 7;,)! - 262TV7;, 

for any < 6 < 1. Taking the supremum over h yields ()4.23|) . ■ 

We now apply this lemma to 1)4. 2U() . with h = /3(— A — /i). Note that 
f{h) = -\j3L^P^{l3,^x) in this case. Let 70 = 7?^ = (1 + z'^ exp(-/3A))-i 
be the minimizer of ()4.21|) . We can infer from (|4.22|) and (|4.20p that 

^(iXin(Xi)Tr(7f -70)' < \PL^ (Po^(/?, ^) - P^[r]) . (4.29) 

If L is large, o^qq small, and z = e^'^ is bounded away from zero, the lower 
bound to the pressure derived in the previous section shows that we can 
restrict our attention to density matrices T with V^\r] > Pq{(3,h) — CagQ, 
for some constant C > 27r. Hence, for such a T, 



(4.28) 



:dX^ n{X^) Tr(7f ' - 70)' < kC(3L^agl . 
Using (|4.23j) instead of (|4.22|) . we obtain in the same way 



(4.30) 



^dX^n{X^ 





Tr(7f -70) 


2 


Tr(7f -70) 


+ Tr7o 



< 2C(3L^aQl 



(4.31) 



By using convexity of the map a; 1— > x^/(x + 1), as well as the fact that, by 
Lemma Tr7o < ^L^qq, (|4.3H) implies that 

^dX^ n{X^) |Tr(7f' - 70) | < L^^QoVCagof^ (l + y^iC^) . (4.32) 
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We thus have an upper bound on both the average Hilbert-Schmidt norm 
of the difference of 71*""^ and 70 and the average difference of their trace. We 
can thus obtain a bound on the average trace norm of their difference, which 
wiU be needed in the next subsection. 

Let a and h be two positive trace class operators, let P be a projection 
with finite rank, and set Q = 1 — P. With || -11^ = (Tr[| • denoting 
the Schatten p-norm, we have 

||a-6||i < ||(a-6)P||i + ||aQ||i + ||6g||i 

< ||P||2||a-6||2 + ||a||f IIQaQllf + ||Q6Q||f . (4.33) 

The trace norm of a can be estimated by II a||i < ||6||i + |Tr (a— 6)|. Moreover, 

IIQaQlli = TTaQ = TT[hQ + {a-h) + {a-b)P] 

< ||Q6Q||i + |Tr(a-6)| + ||P||2||a-6||2. (4.34) 

In conclusion, we thus obtain that 

\\a-b\\i < ||P||2||a-fe||2 + 2(||6||i + |Tr(a-6)|)^/^ 

X (IIQ6QII1 + |Tr(a - 6)1 + ||P||2||a - 6II2)'/' . 

(4.35) 

We apply this inequality, with a = j^'' and 6 = 70, using the estimates 
(|4.3()j) and (|4.32|) . We choose P to be the projection onto the subspace of 
Til = L2(Ai;C) where -A < Kgl^^ for some K > 0. Using Lemma 0J we 
have ||5||i = Tr7o < ^L^qq. Moreover, again by Lemma 01 we can estimate 

\\P\\l = TrP< ^ J dpeiKeT-p") = L'qo-^ (4.36) 
and also 

IIQ6QII1 = TrQb<zTTe^^e{-A-Kgl^^) 

, 1 2/3NT. i/3A . .3 ^exp(- 

< zexp ( - ^(3Ke^ )Tre2'^ < L eo ^2^/3^/^80 

(4.37) 

Similarly to the discussion of the term (|3.45l) in the calculation of the lower 
bound, this last fraction is exponentially small in a^go if we choose K = 
{a'^go)~'^ for some u > 0, uniformly in z for bounded 1/z. By using (|4.3()|) . 
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H4.32() and (|4.35|) - H4.37() . as well as the Schwarz inequality for the integration 
over , we infer that, for small a? q and z bounded away from zero, 



^dX^n{X^ 



7t 



70 



l/3\ 1/4-3I//8 /a 2/3U/4 



(4.38) 



for some constant depending on v. 



4.3.2 Comparing different boundary conditions 

In the following, it will be convenient to compare 71*"^ not with 70 but 
rather with 7per, which is the minimizer of (|4.21|) with h = /?(— Aper — /^), 
Aper denoting the Laplacian with periodic boundary conditions on the cube 
Al. Note that 7per has a strictly constant density. The density matrices 70 
and 7per agree in the thermodynamic limit, however. This can be seen as 
follows. Since the quadratic form domain of A is included in the quadratic 
form domain of Aper, we can use 70 as a trial density matrix of £h with 
h = Aper — /i). Since the pressure is independent of boundary conditions 
in the thermodynamic limit, 

hm J3(/(/3(-Aper-/x))-/(/3(-A-;.)))=0. (4.39) 

Thus Lemma [71 together with (|4.35|) implies that, for fixed (3 and /i, 

lim ^ II70 -7per|li = 0. (4.40) 

4.3.3 A bound uniform in the fugacity 

The estimate 1)4. 38() is not uniform in the fugacity z = e^^. In fact, Pqq^^ 
grows like ln(2:) for large z. Note that large z corresponds to the low- 
temperature limit where the system approaches its ground state. Hence, 
for large z, we will compare 71*"^ with the Fermi sea corresponding to the 
ground state, namely = 9{fi + A). More precisely, we are going to use 
(|4.38() only in the case when Pqq''^ < {a'^ Qo)~^^^ ■ For the case of larger z, 
where Pqq^^ > {a^ Qo)"^^'^ , we now derive a separate bound. 

We start with the following estimate. Let = 1 — and e(^) = 
Tr (—A — fi)Pfj,- For non-negative numbers r, s > 0, and for any operator 7 
with < 7 < 1, 

Tr(-A-^ + rP/,-sQ^)7 

> Tr (-A - ;u + rP^ - sQ^)0{fi + A - rP^ + sQ^) 

= e{fi-r) + e{fj. + s)-e{fj.) + sTrP^. (4.41) 
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Hence 



Tr(-A-^)7-e(/i) 

> [e(/i - r) - e(^) - r TrP^] + [e{fi + s) - e{p) + s Tr P^] 

+rTrP^(l-7)+sTr(l-P/,)7. (4.42) 

Note that, in the thermodynamic hmit (and for ^ > 0), 



1 



for r > 0. Similarly, for s > 0, 

1 1 / 



(4.44) 



Now if we choose r = 47r^/i -"^/^L ^TrP^(l — 7) and s = 0, ()4.42p implies 
that 

Tr (- A - ^)7 > eifi) + 2^V"^^^^"^ [T> ^/.(l " 7)]^ - o(^^) (4.45) 
for any < 7 < 1. On the other hand, choosing r = and 
47r2 TV (1 - P^)7 



+ 27r-i/3L-i (Tr (1 - P^)7) 



1/3 ' 



(4.46) 



a simple estimate yields 



Tr (- A - ,)7 > e(,) + ^- ^lll^^il^ _ _ .(^3) . 

^%V2 + 27r-V3L-i (TV (1 - P^)7)^/=^ 

(4.47) 

Using (|2.4j) . we have 

Tr (-A - ^)7 - ^S[j] > -^P^{f3/2, /i) + ^Tr (-A - . (4.48) 

We use this estimate, together with ()4.45() . on the last term on the right 
side of 1)4. 20() . with 7 = 'j^'' . We restrict our attention again to states with 
V^lT] > Po{P,n) - Cagl as above. Note that e{ii) = -^L^P^{oo, fi). Let 
AP(^(/?, /x) denote the expression 

APo^(/3, /i) = -2Po^(/3, ^) + P(f (/3/2, ^l) + Po^(oo, /i) > . (4.49) 
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The positivity follows from convexity of Pq{i3, h) in 1//3. From (|4.2U|1 . (|4.48|1 
and (|4.45|) we obtain 

-1^ dX^ n{X^) [tV P, (1 - ^) j ' < ^ {ACagl + AP,^{P, ^)) + o(l) . 



47r2 



(4.50) 

By repeating this argument, this time with (|4.47j) in place of (|4.45|) . and 
using convexity of the map x i— > x'^/{l + x^/^), we also obtain the bound 



^^dXin(Xi)Tr(l-P^)7f 



1 + 27r-V3^~i/2 



^fdXin{Xi)TT{l-P^)j^ 



1/3 



<l^{ACagl + APo^{P,^i))+o{l). (4.51) 



1/2 



27r2 



We claim that 



APoHP,^i)< 



37r2 /32 



(4.52) 



in the thermodynamic limit. To see this, first note that APq{P,^) < 
Po^(/?/2,/i) - Po^(oo,/i). Moreover, 

hm [Po^(/3,^)-Po''(oo,/i)] = — ^ / 



+ 



1 



(27r)3/? 



dp In yl + z e 
dp In ( 1 + ze~ 



(4.53) 



Estimating ln(l + x) < x, we see that the first integral is bounded by 



47r//2 



3 z 



1/2 



2tt /xV2 z-1 , 27r m^^' 



< 



3/3 z - 3 /3 



(4.54) 



In a similar way, the second integral is bounded by 



An z 



1/2 



(4.55) 



Hence we arrive at (|4.52|) . 
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Now assume, as explained above, that (ig^^ > {a^Qo) For small 



a Qo, this means that z has to be large. In this case. 



1/2 



27r2 



{ACag'o + AP^{(3,fi)) < const. gt{ago^') 



1/3N2/3 



(4.56) 



Note that, for P a projection, 

h-^'lli < ||(7-l)^l|i + ll7(l-^)l|i 



Moreover, 



< ||P||2||(7-l)P||2 + ||7l/2||^||^(^_p)||^ 

< ||P||2[Tr (1 - 7)P]'/' + II7I/2II2 [Tr7(l - P)] , 



|7i/2||2 = Tr7 < TrP + Tr7(l - P) . 



(4.57) 



(4.58) 



Note that TrP^ < L^iGTr"^)-^ fi^/^ (using Lemma 0} . By combining the 
estimates KEU^ . (|13T]) and (|T^ - (II3H1) . with P = and 7 = 7:^^ 

obtain that, for small a'^go and Pgfj^^ > {a^ go)"^^'^ , 



we 



L3 



dX^ n{X 



P, - if 



< const. (a^Q 



1/3U/6 



(4.59) 



This inequality supplements 1)4. 38(1 in the case of large z. 

It remains to estimate ||P^ — 7o||i = TrP^(l — 70) + Tr7o(l — P^) 
the thermodynamic limit, 

1 



lim 7^||P/ 



In 



70 111 



< 



< 



(27r)3 
1 

(2^ 

1 /iV2 



dp 
dp 



e{^x-p') ^ 9{p^-fi) 



1 + ze-(^P 1 + z-^e(^P 



1 + 



ze-f^P 
1 



-lg/3p2 



2/3// 



(4.60) 



67r2 f3 

where the last inequality is derived in the same way as in (|4.53|) - (|4.55|) . In 
case Pgf)^^ > {a^ go)~^^^ , as considered here, this last expression is actually 
bounded by (agQ^^)^^^ , and therefore negligible compared with the error 
term on the right side of (|4.59j) for small a'^go- 

4.3.4 Summary of this subsection 

To summarize, we have shown in this subsection that for a density matrix 
r satisfying P^[r] > Pq{P,^) — Cag^ we have the a priori bound 



dX^ n{X 



7per - 7t 



< CuL^go{ago 



+ o{L^) (4.61) 
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for > 0, for some constant C^, depending only on v (and C above), but 
not on z. Here, 7per denotes the one-particle density matrix of a system of 
non-interacting fermions (at inverse temperature (3 and chemical potential 
/i) with periodic boundary conditions on the cube A.^. 

As a side note, we remark that (|4.61|) implies, in particular, that the 
reduced one-particle density matrix of the dilute interacting system, 7^, is 
close to the one for non-interacting particles. This, in turn, proves the 
inequality H1.8() in Corollary^ with a worse error term than the one given 
in however. 



4.4 Putting Things Together 



We now show how the estimates of the preceding subsections can be com- 
bined to prove the desired lower bound on the pressure, given in Theorem^ 
Let be the one-particle operator given in 1)4. 1U() . It follows from (|4.11() . 
1)4. 17() and 1)4. 18() that, for any density matrix F, 



1 



+{l - 6){l - >c)i:dX^ n{X^)TTn, Wxilf .(4.62) 



Here, we used that T is symmetric with respect to exchange of the | and 
j-particles, by assumption, which implies in particular that 7^ = 7|. With 
the notation introduced in Lemma [3 



Tvn, W - 7t] - ^S[7t] > -^f{m - ^)) 



In the last term in ()4.62|) . we use that 



Tr 



Hi 



Wxnf^ > Tr^i [Wxnper] - WWxlWochf' -7per||l 



(4.63) 



(4.64) 



We have already estimated || Wj^i ||oo in 1)4. 16() . and this estimate is indepen- 
dent of XK Moreover, the X^-average of ||7|'"^ — 7per||i bas been estimated 
in (|4.()H) . Because of translation invariance, 7per has a strictly constant 
density, given by L-^Tr 7per. Thus 



L3 



[0,L]- 



dx Wxi ix) . 



(4.65) 



We note that L ^Trjper = ^go + o(l) in the thermodynamic limit. To 
estimate the integral in (|4.65|) . we write Wxi = W^ — W^, where W_|_ 
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denotes the positive terms in H4.5() containing f7, and —W- the negative 
ones containing wr. For the negative part, we can use (|4.13j) to get the 
upper bound 

f aE? 

/ dxW-{x) < const. — (4.66) 
J[0,LP 

Here, \X^\ = M denotes the number of spin down particles. Note that 

^dX^n{X^)\X^\=TT-ii. (4.67) 

This trace can be estimated using (|4.61|) . which in particular implies that 

iL-^TrT^^ _ < CL^Qo^agy^f''''' + o(L=^) . (4.68) 
Moreover, using J U = 4Tr, 

/ dxW+{x) = (1 - e)a / dxU{x — xl.) 



{k-.xi&x^^} 
> (1 - e)4-iTa 



Xj^l — const. -;::7T 



(4.69) 



for any fixed XK Here, \xj^\ denotes the number of elements in X^, i.e., 
the number of x^'s in X^ whose distance to the nearest neighbor among the 
x^'s for / 7^ A; is bigger than 2R. The last term in square brackets bounds 
the number of x^'s that are closer than a distance R to the boundary of 
the box. Since the distance between the x^'s in X^ is bigger than 2R, the 
number of such x|,'s close to the boundary is bounded by const. !?■ /R^. 
We now need an estimate on Note that 

1 

\x\,\>\X^\-{2RfY.^^ (4.70) 

k=l k 

where denotes the distance of xj, to its nearest neighbor in X^. The last 
expression can be bounded from below using the operator inequality 

M M 
k=l k 
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which holds on the anti-symmetric tensor product /\^^ L^(M'^), and is proved 
in [HI Thm. 5]. Here, c is some positive constant satisfying c < 48. We thus 
have that 

^dX^n{X^)\X^^\>TTn,{l + c{2RfA)ji. (4.72) 

The negative last term involving the kinetic energy can be canceled by an 
appropriate choice of 6 in (|4.62)) . 

By choosing 6 = 2'iracgo{2R)'^ and taking the thermodynamic limit L 
oo in (|4.62|) . the above estimates imply that 

P{P, fi)<- hm sup -^f{P{h^ - ^)) - 27Tael {I - 5 - e - k) 



+CuaQl 



(4.73) 



It remains to estimate the first term on the right side of this expression. 
To this end, let T(p) = (1 - 5)p^{l - (1 - >i)x{pf) (compare with (BTTnll V 
We claim that 

hm^inf -^f{m - m)) > dp In (l + zexp ( - /3T(p))) . (4.74) 

This can be seen using coherent states as follows. Let n be a real function 
supported in a cube of side length I, with / dxu{x)'^ = 1, and let Cy,k{x) = 
u{x — y)e*^^. Denoting by \iy^k){iy,k\ the projection onto ^y^k, we then have 
the resolution of identity (2t:)~^ J dkdy\^y^k){Cy,k\ = 1- Since x i— > — xlnx 
is concave for < x < 1, this implies that, for any fermionic one-particle 
density matrix 7, 



Sh] < -{2Tr)-^ j dk dy [0{y, k) In Q{y, k) + {l - 0{y, k)) In (l - g^y, k))] , 

(4.75) 

where Q{y,k) = {6.y,k\l\Cy,k) ■ Note that < Q{y,k) < 1, and that Tr7 = 
(27r)^3J dkdy Q{y,k). With 7 (A;) = (27r)^3J dx' 7(x, x') exp(iA;(x' — x)) 
denoting the Fourier transform of 7, we can write 

Tih'^^j = j dpT{p)^{p) . (4.76) 

Moreover, a simple calculation shows that 

/ dk dy T{k)giy, k) = {2-nf j dk dq T{k + q)7ik)\u{q)\'^ . (4.77) 
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Note that j{k) > 0. In the case considered here, T has a bounded Hessian, 
i.e., didjT{k) < 2C as a matrix, for some constant C (depending only on 
the parameter s used in the definition 1)4. 12() of x)- Therefore T(k + q) < 
T{k) + qVT{k) + Cq^. Inserting this estimate into (|4.77|) and noting that 
/ dqq\u{q)\^ = since u is assumed to be real, we obtain 

Tr [/i^7] > (2^)-3 / dk dy g{y, k) ( T(A:) - C [ dx | Vn(x)n . (4.78) 



Now choose u such that J dx\Vu{x)\'^ < C'i ^ for some constant C inde- 
pendent of L Then (|I?7^ and (|I?75|) imply that 

Tr[(/iX_^)^]_l5[7] 



> (27r)-3 j dk dy {T{k) - - CCT'^) Q{y, k) 

+ (2^)3^ j '^^ [^^y^ ^) Siy, k) + {l- Q{y, k)) In (1 - Q{y, k))] . 

(4.79) 

Note that g(y, k) ^ Q only if y is an cube of side length If we minimize 

the integrand for each fixed y (in this cube) and fixed /c, we see that the 
right side of (|4.79j) is bounded from below by 

_(L + 2^)3-J_ [ dk\n(l + ze^cc'/i\-prik)\ (4 

{2TTy[3 J \ ) 

Dividing by L^, letting L — > 00 and then £ — > 00 we arrive at (|4.74|) . 

It remains to estimate the integral on the right side of (|4.74j) . and com- 
pare it with — iPo(/?)/")- We do this in two steps. First, note that 



2 



(2^)3/3 



y dpln(l + ze-'3T^(^')) -(l-5)Po(/3(l-<5),/x/(l-<5)) 



2 1' 1-1- z(=~^'^(P) 

^dp\n ^ ,t ' -.w (4-81) 



(27r)3/3 7 ^ 1 + ^e-Z^Ci--?)?' 

By the definition ()4.12|) of X) T(p) = (1 — 5)p'^ for |p| < 1/s. Hence the 
integrand in (|4.81j) is only non-zero for \p\ > in which case T{p) > 
(1 — 6)xp'^. Hence (|4.81() is bounded from above by 

dp In— -57T-7^^^ < TTT-TT^S / dpe '>^P 



(2^)3/3 - 1 + ze-/3(i-%^ - (2^)3/3 

- y2vr2/35/2(l-25)x' ' ^^"^^^ 
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where we estimated \p\ < (2/3) ^/^exp(|/3p^) in the integrand in order to 
evaluate the last integral. Secondly, by simple scaling, 

(1 - 6)Po (/3(1 - 5),fi/il - 6)) = (1 - 5y^/^PoiP, fi) . (4.83) 

We thus arrive at the following estimate: 

P(/3, ^) < (1 - dy'/^PoiP, /i) - 27ragl {I - S - e - x) 



^ , /, ,1 J 1 1 n /' l/3U/6-i. 



-,+ ^l + -max|^^3-^,-^|)(a,o 



V27r2^5/2(i_25)x' ' ^^-^^^ 

with 5 given as above hy 6 = 27racQo{2R)'^ . It remains to choose s, R, e and 
X. We take 

D -1/3/ 1/3n1/22 -1/3/ l/3\l/66 / 1/3n1/33 ^r-x 

and 

^=^aQl^Y'''"'- (4.86) 
Note that with this choice of parameters, 5 = 8TTc{a^ qq)^^^^^ . Hence 

6Po{P,fi) = const, aeli^' go f ^ff^^L ^ (^.87) 

and the last fraction is uniformly bounded for z bounded away from zero 
(as already argued in the lower bound after (|3.47() ). Moreover, since xs~^ = 
{o^ Qq)~^^'^ , we see by the same reasoning as in the lower bound in ()3.45p 

1 /3 

that the last term in (|4.84|) is actually exponentially small in ag^ , uniformly 
in z for bounded \/z. Inserting ()4.85p and ()4.86|) into 1)4.84(1 above, we thus 
obtain 

P(/3, /i) < Po(/3, - 2T:agl (l - C,{z) {agl''') (4.88) 

for any > 0, with Ci,{z) bounded for bounded 1/z. This finishes the proof 
of the upper bound in Theorem ^ 

5 Proof of Corollaries [T] and [2] 

In this final section, we show how to derive Corollaries Q and |2 from The- 
orem n We start with proving Corollary ^ The essential ingredient is 
convexity of P{f3, fj.) in /i. It implies that 

/. /. P(Aj^ + 5)-PiM 
g-{P,fM) < g+{P,ti) < (5.1) 
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for any 6 > 0. Using (|1.7() as well as the fact that QoiP,^) is monotone 
increasing in this yields 

Q+{P,H) < ^ (5.2) 

+]qo{P, + Sf/^C^e^^^) {agoiP, /x + 5)^3)1+" . 
Again by convexity, Po{P, fj, + 5) — Po{(3, /i) < Sgo{P, ^ + 5). We choose 

5 = {aQ^{f3,iJi)^'^f'"'^'^ max{^, . (5.3) 

With this choice of 5, it is then not difficult to see that 

Qo{l3,ll + 5) , l/3N(l+a)/2 ,^ .s 
< 1 + const, [ag^ y (5.4 

for some constant independent of ^ and (3. Using this estimate, ()5.2() implies 
that 

Q+{(3,^x) < 0o{P,^^) (l + C«(z)(a£,y=^)(^+")/') (5.5) 

for some Ca{z) that is uniformly bounded for bounded 1/z. 

A lower bound on f3_(/3, /i) can be obtained similarly, using that, for 
6>0, 

PiP^j^-6)^P(M .... 
g-iP,fj') > ^ • (5.6) 

Proceeding along the same lines as in (|5.1j) - (|5.5j) above, this proves Corol- 
lary [H 

Next we prove Corollary [JJ For ^ > 0, let /(/3, g) = sup^[fig — P{j3,ii)] 
denote the free energy per unit volume. If the supremum is attained at some 
fi, then convexity of P{f3, fi) in fj, implies that 

g-{(3,fi)<g<g+{p,fi). (5.7) 

Let /o(/3, g) = sup ^[fj, g— Po{(3, /x)] denote the free energy density for the ideal 
Fermi gas. The supremum is achieved at fiQ, determined by £>o(/?)/"o) = Q- 
Hence we immediately get the lower bound 

/(/?, e) > f^oQ - PW, /io) > UP, g) + 27rag^ - {e''^')ag^ {ag^^T ' (5-8) 

where we used (|1.7() to estimate P(/3, /xq) in terms of Po(/?, /xq)- 

To get an upper bound on the free energy, we first make an a priori 
estimate to ensure that /x is close to ^o- Suppose that /x < /xq. Then, by 
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H5.7|) and monotonicity of Q+{l3,fj,) in fi, g = Qo{(3.,iio) < q+{I3,iiq — 5) for 
all 5 < ^0 — A- Using (|1.8|) . this implies that 

eo(/3, A'o) < f?o(/?, /^o - 5) (1 + /^o - 5)) (5.9) 

for all (5 < ^0 ~ A) where we denoted 

=a„(e^'^)(a£.o(/3,/x)'/=')^'''°^^'. (5.10) 
Ineq. 1)5. 9(1 then implies that 

5 < const, (ae^/^)^^^"^/^ max{^o, (5.11) 

for some constant independent of g and [5 (compare with (|5.3|) ). Denoting 
the right side of ()5.11|) by 5, we therefore see that fl > fJ-o — 6. 

The same method works in the case when fi > fiQ. One uses that 
go{P,Ho) ^ Q-{f3,l^o + ^) for all 5 < fi — fiQ. Proceeding along the same 
lines, this implies that fi < fiQ + 5. In particular, 

f{/3,g)= sup • (5.12) 

Using (|1.7j) and (|1.8|) . it is not difficult to see that 

P(/3,/i) >Po(/3,/^)-27ra£^o(/3,M)'(l + C„(e^^o)(a£)o(/3,/^)'/')") (5.13) 

if l/U — //o| ^ for some Ca(z) that is uniformly bounded for bounded 1/z. 
Inserting the bound 1)5. 13(1 into (|5.12p proves the desired upper bound on 

f{P,g)- 

Note that e^^° is bounded away from zero for bounded l/{Pg'^^'^). Hence 
(|5.13|) and (|5.8j) imply the statement in Corollary 121 
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